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This report outlines the theory used in FASTER -III, a Monte Carlo 
computer program for the transport of neutrons and gamma rays in com- 
plex geometries. The code has the additional capability of calculating 
the minimum weight layered unit shield configuration which will meet a 
specified dose rate constraint. It includes the treatment of geometric 
regions bounded by quadratic and quadric surfaces with multiple radia- 
tion sources which have a specified space, angle, and energy dependence. 
The program calculates, using importance sampling, the resulting number 
and energy fluxes at specified point, surface, and volume detectors. 

Results are presented for sample problems involving primary neu- 
tron and both primary and secondary photon transport in a spherical 
reactor-shield configuration. These results include the optimization of 
the shield configuration. 


Section 1 

INTRODUCTION AND SUMMARY 

The original FASTER program (ref. l) contained 
a number of new techniques which provided the capa- 
bility of obtaining accurate radiation levels at 
specified points in complex geometries. Prior use 
of FASTER indicated a need to broaden the overall 
program capabilities, automate the importance 
sampling, increase the computational efficiency, 
and revise the users manual. This revised program 
has been designated FASTER-III to distinguish it 
from earlier versions. 

A specific program capability permitting the 
calculation of minimum weight layered unit shield 
configurations for mobile nuclear reactor applica- 
tions, e.g. , nuclear propulsion for aircraft, sur- 
face effect vehicles, and spacecraft has recently 
been developed. The basic Monte Carlo transport 
method was extended to include a calculation of 
partial derivatives of the radiation fluxes with 
respect to specified shield dimensions. These de- 
rivatives are then used to define exponential re- 
lationships used in the shield optimization pro- 
cedure. This optional program feature is described 
more completely in Section 2. 

Data preparation is simple, with very little 
judgment required to set up the importance sampling 
for most problems. The code also has a unit shield 
weight optimization capability. 

Particularly noteworthy features of FASTER- 
III are the following: 

(1) A calculation of optimal importance sam- 
pling parameters based on partial derivatives of 
the variance (Section 2.3). 

(2) The acceptance of data in either fixed or 
variable field formats including the ANISN-DTF for- 
mat for neutron cross sections. 

(3) The calculation of time-dependent neutron 
and photon transport (using time moments and/or 
time intervals) including an optional exponential 
atmosphere. 

(4) The improvement and addition of importance 
sampling models with the various importance sam- 
pling parameters built into the program. 

Various program features are described in 
Refs. 2 to 6. 

The application of the FASTER-III program to a 


shield optimization problem is discussed in Section 
3. The problem involved a spherical re actor -shield 
configuration and included primary neutrons and 
both primary and secondary photons. Conclusions 
and recommendations are presented in Section 4. 

Section 2 

ANALYSIS 

The techniques used in calculating optimum 
shield configurations and optimum importance sam- 
pling parameters are summarized below. The dis- 
cussion is given in three parts: dose rate deriv- 

atives with respect to shield layer thicknesses, 
optimization procedures, importance parameter op- 
timization. 


2.1 Dose Rate Derivatives 

The dose rate at a point detector for a 
specified reactor shield configuration is written 
as : 


J 

D(z) = KjCpjCi) (1) 

J=i 

where J is the total number of energy groups for 
both neutrons and photons (including secondaries), 
cPiCx) is the particle flux in the jth energy group, 
and Ri is the response function to convert from 
flux to dose rate. The rate of change of the dose 
rate with respect to a shield layer thickness is 
simply 
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where L is the total number of shield layers and 
t 1 is the thickness of the 1th layer. The equa- 
tion used by the program for determining the flux 
is written as: 
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where N is the total number of histories tracked 
via the Monte Carlo method, k is the number of 
particle collisions, is the position of the 

kth collision of the n th history, Sjj^u^) the 
number of particles in the jth energy group 
emerging from in the direction uj^ of the 

detector per unit solid angle, and KjX^kn^ rep- 
resents the material and geometric attenuation 
kernel for particles in the jth energy group 
going from to the detector. 

The partial derivative of the flux with re- 
Zth shield layer thickness is 
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The summations are a minor part of the calculation. 
Therefore, the notation is simplified by concentra- 
ting on the elements in the summation 



where 0jkn represents the contribution to the 
flux in the jth energy group from the kth col- 
lision of the n th history. This equation is re- 
written as 
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The partial derivative of the partial path 
length Sjjj with respect to the shield layer thick- 
ness tj is zero unless the mth region tra- 
versed is affected by a change in t In partic- 
ular, if t^ is a characteristic dimension of the 
region, i.e., its thickness, then 

= ~ '> »*knm = Hkn * 2lmm ( 10 ) 
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where ji knin is the cosine of the angle measured 
from the surface normal r^mmy w i*th which the par- 
ticle crosses the boundary of the region. 

In the strict sense, the change of the thick- 
ness of one shield region can affect other shield 
regions. In particular, for a spherically symmet- 
ric reactor-shield configuration, an increase in 
the thickness of a shield region forces a movement 
of all shield regions having a larger radius. The 
inclusion of these effects in the above equation 
unnecessarily complicates the analysis and the cal- 
culations. The primary effect of changing a shield 
region dimension is to change the number of mean 
free paths which particles have to traverse in 
reaching the detector. Therefore, in calculating 
the derivatives, only the effect of the material 
attenuation is treated. 

The derivatives at a specific boundary cross- 
ing m' then simplify to: 
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where M is the total number of regions traversed 
from z^ to the detector, s is the path length 
for the mth region traversed, is the total 

cross section of this region for particles in the 
j th energy group, and s is the total distance 
from z^ tp the detector, i.e., 

M 

s - E s ” 

m=l 

A substitution of this kernel gives: 


= - a jm'Aw 

where m' is the index of a region having tj as 
a dimension. The partial derivatives of the par- 
ticle weight with respect to the shield dimensions - 
the first term in brackets in Eq. (6) - are zero 
at the point of origin of all primary particles. 

For subsequent particle collisions, the deriva- 
tives are calculated using the relationship be- 
tween particle weights on subsequent collisions: 
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where S* ^y^) is the number of particles 
coming out of ' * the previous collision point in the 
direction y^ n and in the ith energy group, 
^i(5k-l n> ^kn) the attenuation kernel between 
particle collision points, TijC^kn, v^ n • U}^) is 
the scattering kernel for transfer of particles from 
group i to group j, and p£n(zkn) is the Probability 
density function used in selecting the collision 
point. 

A straightforward substitution gives 
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where lJ n (Zkn^ ts a probability density function 
used to select the particle direction, s = \hrn ~ 

5k -1 nl is the distance of the selected collision 
point from the previous collision point, A(s) is an 
importance factor for each region which changes 
discontinuously at region boundaries, and a(s) is 
an effective cross section which changes discontin- 
uously at region boundaries and which may change 
continuously within a region. 


The derivative of the logarithm of pj^z^) 
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After some manipulation, this reduces to 


involves only those terms which change when a 
shield dimension changes, i.e., 
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Let Sj denote the distance to a boundary in- 
volving the Ith shield dimension. If the first 
term on the left side of Eq. (17) is affected by a 
change in this shield dimension, i.e. if s > s 7 , 
then 
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The first term in brackets in Eq. (14) is the 
same partial derivative for collision k-1 as the 
partial derivative now being calculated for colli- 
sion k. Therefore, it is known, either identi- 
cally zero for k=0, or as determined from Eq. (14) 
for k > 0. The second term in brackets in Eq. (14) 
is similar to the second term in brackets in Eq. (6) 
and is therefore determined by Eq. (ll). The last 
term in brackets involves the definition of the 
probability density function used to select the 
collision point z^ur 

The probability density function for a colli- 
sion point has the form 


where a(s^) is the effective cross section at the 
boundary of the shield and p^kn i s the cosine the 
particle path makes with the outer shield normal. 

If there is any crossing involving the 1th shield 
dimension, the second term in Eq. (18) will always 
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Curved shield surfaces may "be crossed more than once 
along the path between two particle collision 
points. Therefore, a summation of Eqs. (18) and 
(19) over every intersection involving the 7th 
shield dimension is required to completely evaluate 
Eq. (17). 


2.2 Optimization Procedures 


The shield optimization calculation yields the 
set of shield layer thicknesses t' = (t-j, t£, . .., 
t|, ... t£) such that the dose rate, D(t*), meets 
the dose constraint. The Monte Carlo calculation is 
performed for an initial set of shield layer thick- 
nesses t = (t^, ... t^, . .., t^) and yields a 

set of fluxes, cpj(t), j = 1, 9 , ..., J and deriva- 
tlves, dcpj(t)/dt z , j = 1, 2, •••, J; 1 = 1, ?, •••, 
L. The assumption is made that the fluxes vary 
exponentially with respect to shield dimension 
changes in the form 

cpj(t') = cpj(t) expQtj • (t 1 - t)] (20) 

where a^ = (a^, aj 2 j •••> follows that 
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In particular 


ckpi(t) 

= aj jcp.(t) (22) 

or 
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( 23 > 

The weight is also expressed as a function of 
the shield layer thicknesses. The weight is de- 
noted hy W(t’) and for spherically symmetric 
shields : 



where is the density of the 7th shield re- 

gion ana r 0 is the minimum shield radius. 

The purpose of the optimization procedure is 
to minimize the weight W(t ' ) subject to the dose 
rate constraint D(t') = D 0 where D 0 is a speci- 
fied dose rate. At this optimum, a small weight 
perturbation in any layer causes the same dose rate 
change. The rate at which dose rate changes with 
respect to a shield weight change in the 2th layer 
is given by 
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The necessary derivatives are: 
dD ( t * ) = 


= ^ Rjaj2tpj(t)exp[aj • (t' - t)J (26) 
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and for spherically symmetric shield: 



(27) 


In arriving at the optimum shield, the total 
shield weight is built up in increments of weight 
£W. Each increment in shield weight is always 
associated with a particular shield layer thickness. 
At each iteration, the particular shield dimension 
is selected by examining the values of the shield 
weight quality factors, Q,. Each factor rep- 

resents the approximate change in dose rate per 
unit change in weight corresponding to a change in 
the 7th shield dimension. Negative ’ s are the 
most usual and correspond to shields for which an 
increase in weight - and shield dimensions - gives 
a decrease in dose rate. Positive * s can occur, 
however, and correspond to shields for which an in- 
crease in weight also increases the dose rate. 

If, at a particular iteration, the dose rate 
is above the dose rate constraint, the minimum 
shield weight increment would correspond to the 
least positive value of those ' s for which 

> 0 and for which tj > t^(min), where t^(min) 
is the minimum value of the 2th shield layer 
thickness. If such a Hi exists, the dose rate 
can be decreased while also decreasing the shield 
weight the maximum amount. If there isn't such a 
Q^, the next best procedure is to find the most 
negative of the Q? ' s for which < 0 and for 
which t| < t^CmaxJ, where t^(max) is the maximum 
value of the 7th shield layer thickness. A 
change in that would give the maximum decrease 

in dose rate per unit increase in weight. 

If the dose rate is below the specified dose 
rate at a particular iteration, the minimum shield 
weight increment would correspond to the least 
negative of those ' s for which < 0 and for 
which tj > t^(min). If such a exists, the 

dose rate can be increased while decreasing the 
shield weight the maximum amount. If there isn't 
such a Q|, the next best procedure is to find the 
most positive of those ^ ' s for which Qj > 0 
and for which t j < t j (max) . A change in that 
would give the maximum increase in dose rate per 
unit increase in weight. 
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Assuming a particular value of the ’ s 

is selected through the above argument s, the cor- 
responding shield dimension t^ is changed by a 
maximum amount At where At is calculated as 
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If this change would put tjJj outside one of 
its specified limits, the value of t^ would be 
set to that limit, i.e., t^min) < t^ ^ t^max) . 
The shield weight increment AW is calculated as 
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subject to the constraint that |AW| < AW Q where 
A^ 0 is a specified maximum shield weight incre- 
ment per iteration. Note that ZW, and therefore 
At , may be positive or negative depending on the 
value of and whether the dose rate is above 

or below the dose rate constraint. 

Once a shield layer thickness is changed, the 
dose, weight, and their derivatives are re- 
evaluated and the entire process is repeated. The 
optimization would be discontinued in several 
ways. If the dose rate equals the dose rate con- 
straint within the relative error of the original 
Monte Carlo dose rate calculation, the program will 
proceed to the next problem - which may be identi- 
cal except with more histories to tighten the con- 
vergence of Monte Carlo calculations. Similarly, 
if all shield layer thicknesses have reached their 
minimum or maximum values, and if the optimum 
shield cannot be determined with these constraints, 
the program would again proceed to the next prob- 
lem. Finally, if the dose rate and dose rate con- 
straint are decades apart in value, the program 
would reevaluate the fluxes and their derivatives 
by Monte Carle every time the dose rate changed by 
more than a specified factor during the optimiza- 
tion procedure. 


Thus the partial derivative of the relative 
error with respect to the parameter a is: 
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The dose rate from the n th history is given by 
J 


D n = E E j£ *3*a. 
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where J is the total number of energy groups, k 
is the number of particle collisions, Rj is the 
flux to dose rate conversion factor for the jth 
energy group, and cpj kn is the flux in the j th 
group from the kth 0 colli si on of the n th his- 
tory. Since 
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2. 3 Importance Parameter Optimization 


The optimization of the importance sampling 
must be performed for some function, e.g., dose 
rate, of the energy-dependent fluxes since there 
is a different optimum for every initial particle 
energy. Therefore, assume that a minimum variance 
calculation of the dose rate is required where 

= N S ° n (30) 
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where N is the total number of histories and D n 
is the dose rate from the n th history and Dpj 
is the average value of the dose rate after N 
histories. The relative error of this dose rate is 
given by 
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Taking the logarithm of this equation and 
then performing a formal calculation of the partial 
derivative with respect to an unspecified param- 
eter a yields 


the calculations required to evaluate Eq. (33) all 
involve the summation of terms which involve 


dD„ 
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The remainder of the analysis, therefore, can 
be concentrated on the partial derivatives of the 
fluxes. All other operations which must be per- 
formed are given above. 

The fluxes typically depend on the detector 
position y, so the equation for the particle flux 
is written as 

<Pjkn(z) = S jkn(Hkn) K l5kn> Z) ( 37 ) 


The transport kernel Kj(z^ n , y) does not in- 
volve any importance sampling parameters so that 


d^jkn^ _ d S jkn(-kn) K j 0^* l) 
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This equation can also be written as 

= S jkn(Hkn)^j(^kn'£) 7“ ^jkn^Hkn) (39) 
da da 

Without going into great detail, it turns out 
that the particle weight S^ kn (u kn ) is composed of 
a purely analytical numerator, v jkn(Hkn) 811(1 a de “ 
nominator which is the product or all the probabil- 
ity density functions used to select the collision 
points, i.e., 
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Since V'jkn(ukn) does not explicitly involve 
any importance parameters, it follows that 


to a set of parameters which can have a reasonably 
simple role. These parameters consist of the rel- 
ative importance I r of each region. Normally 
these parameters are all equal. However, in asym- 
metric problems, it turns out that some regions are 
much more important in terms of their scattering 
contributions to a detector. Therefore, these im- 
portant regions have a larger value of I . 

The region importance enters into tEe selec- 
tion of a collision point through the following 
probability density function: 
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where r is the region in which the collision 
occurs (selected at random), p*(s) is the piecewise 
continuous probability density function in this re- 
gion at the selected collision point (a distance s 
from the previous collision point), H is the total 
number of regions in which the collision could have 
occurred, and is the integral of p^(s’) over 

the partial path length in region h. 

Calculating the logarithm of each side of the 
equation yields: 
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The partial derivative of Eq. (47) with re- 
spect to the specific importance parameter Ig - 
the relative importance of region g - yields^ 


Therefore, Eq. (39) can be re-written as 
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Moreover, the partial derivatives are energy- 
independent so that Eq. (36) becomes 
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The evaluation of the partial derivatives of 
the probability density functions can be written 
as 
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where =0 if region h is not region g and 

5 gg = 1. 

Thus Eq. (48) is evaluated during the random 
selection of the k th collision point and the 
final term necessary to evaluate Eq. (45) and all 
preceding equations has been determined. 

The above analysis is used to calculate the 
partial derivatives of the relative error of the 
dose rate with respect to the relative importance 
I r of each geometric region, and a similar analy- 
sis is performed for the other importance sampling 
parameters. The result of the complete Monte Carlo 
calculation is a set of partial derivatives which, 
for the region importance, are given by 


+ ^ lnp kn ( 5kP (45) 

At the k th collision, the first term on the 
left side of Eq. (45) is known, identically zero 
if k = 0. Therefore, the analysis is completed 
after examining the calculation of the second 
term. 

At this point it is necessary to identify the 
particular importance parameter a. Since most of 
the importance sampling parameters have fairly 
involved roles, the technique will be applied here 



(49) 


where bD n /bl r is obtained from Eq. (44) using 
Eqs. (45) and (48). 

After the calculation is completed, optimal 
values of the importance sampling parameters are 
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calculated by requiring that the relative error be 
zero - not actually achieved of course. 

By a first order expansion 



z=l 


where R is the total number of regions. A simple 
gradient analysis says that 1^ - I z should be 
proportional to &E N /c)I z so that “ 
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where , by substitution into Eq. (50), 
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where z is a general phase 
space coordinate 
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The program prints the optimum values of I z 
and other importance parameters after completing - 
the Monte Carlo flux calculation. This analysis 
is performed for every response function. After 
more experience is obtained with the technique, the 
program could be modified to change these param- 
eters internally corresponding to a specified re- 
sponse function. 


Section 3 

SAMPLE HI OB LEM RESULTS 


Two problems were investigated using the 
shield optimization capabilities of the FASTER -III 
program. Both problems involved a spherical 
reactor-shield configuration and included primary 
neutrons and both primary and secondary photons. 

The two problems were similar except for the 
power level, 375 MW and 600 MW respectively. Both 
problems used a flat radial distribution for the 
primary neutron and photon source distribution. 

The primary photon source included an infinite 
operation equilibrium fission product term. 

The core radii for the two problems were 82.38 
and 96.38 cm respectively, corresponding to a power 
density of 4.53 MW/ft 3 . Following the core was a 
7. 62 cm Be reflector; a 5 cm depleted uranium 
shield; three depleted uranium-borated water shield 
layers of 57, 15, and 15 cm thickness and 6.4, 

4.6, and 2.8 gm/cm 3 density respectively; and a 
117 cm borated water shield. This base line shield 
configuration was based on parameters obtained from 
SANE-SAGE calculations and subsequent calculations 
using the UNAMIT program, Ref. 7. The reactor- 
shield compositions are given in Table 1. 

The primary neutron transport calculation 
utilized multigroup cross sections for 26 energy 
groups. Fifteen energy groups were utilized for 
both primary and secondary photons. The secondary 
production cross sections included both inelastic 
and capture gammas. 

These initial configurations were each analy- 
zed for a point detector 30 ft from the core cen- 
ter by following approximately 500 energy- 
dependent packets of primary neutrons and photons 
and approximately 7000 packets of secondary pho- 
tons. The dose rates obtained from these calcu- 


lations are tabulated in Table 2 including a break- 
down by secondary source region. Each of these 
problems required about 28 minutes on the UNIVAC 
1108 computer. 

The basic calculated dose rates and dose rate 
derivatives were also used by the FASTER -III pro- 
gram to calculate the minimum weight shield con- 
figuration which would give a dose rate of 0.25 mr / 
hr at the specified detector point. The final 
shield configurations following the optimization 
are given in Table 3. 

In both cases, the optimum shield configura- 
tion is significantly different from the base line 
configuration. Since the base line configuration 
was not generated by the FASTER -III program it is 
difficult to discuss many factors entering into 
that calculation which would account for the dif- 
ferent optimal configuration. It is noted, how- 
ever, that the base line configuration was gener- 
ated using parameters corresponding to a calculated 
dose rate an order of magnitude below the specified 
dose rate constraint, Ref. 8. As such, the base 
line configuration used in the FASTER-III program 
was determined from an extrapolation of a different 
base line configuration. 

A more critical critique can be made of the 
FASTER-III results independently. First it is 
noted that neither problem saw a significant con- 
tribution, less than a few percent, from photon 
sources in the core region. In fact, the 600 MW 
reactor dose rate from this source was about a fac- 
tor of two less than it was for the 375 MW reactor. 
This difference is ascribed to the problem statis- 
tics since core photon sources see approximately 
30 mean free paths of shield material. Therefore, 
it is doubtful if this dose rate component is con- 
verged within a factor of two after only 500 pack- 
ets but this does not introduce a significant error 
since the original contribution was only two per- 
cent of the total dose rate. 

The small contribution from core photon 
sources decreases the amount of high Z shields 
required around the core. Therefore, both problems 
gave a significant change in the first two shield 
dimensions during uhe optimization. In the 375 MW 
problem, the first mixture of depleted uranium- 
borated water (p = 6.4 gm/cm 3 ) was eliminated en- 
tirely. In the 600 MW problem, the depleted ura- 
nium and most of the first mixture were eliminated. 

The main difference between the two FASTER-III 
calculations was the shift in the placement of 
lighter shield mixes towards the core for the 
600 MW problem. An examination of the secondary 
photon dose components indicates that the contri- 
bution from the outer two shields was about 25 per- 
cent for the 375 MW reactor and almost 50 percent 
for the 600 MW reactor. Since these sources de- 
pend on the neutron attenuation through the closer 
regions and since lower effective Z materials 
are better neutron attenuators on a weight basis, 
the 600 MW problem tends to replace high effective 
Z material with a lower effective Z material. 

The differences in the contribution from 
secondary sources in the outer shield regions is 
greater than expected for the nominal difference 
in the core region. Therefore, much of the dif- 
ference in these sources must be ascribed to sta- 
tistical variations. In fact, both problems had 
approximately 25 to 30 percent calculated relative 
error in the total photon dose rate. It should be 
noted that the FASTER-III program includes a num- 
ber of importance sampling techniques which could 
be used to decrease this error. However, both 
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problems were run using the "built-in definitions of 
importance parameters. Alternatively, more his- 
tories could have been used although the computer 
time requirements would have become excessive. 

Section 4 

CONCLUSIONS AND RECOMMENDATIONS 

The FASTER -I II program was developed to calcu- 
late neutron and photon fluxes at specified points 
in complex geometries. Alternatively, it can also 
calculate fluxes averaged over specified surfaces 
and volumes. The program was designed such that 
data preparation is simple and so that very little 
judgment is required to set up the importance sam- 
pling for most problems. The FASTER-III program 
satisfies these requirements very well. 

The shield weight optimization capability in- 
cluded in the FASTER-III program permits the cal- 
culation of both base line radiation levels and op- 
timal shield thicknesses all in a single computer 
run. However, the very large attenuation factors 
involved in the demonstration problems yielded some 
questionable results. In particular, the statisti- 


cal differences in the relative contribution from 
various secondary source regions caused correspond- 
ing variations in the relative distributions of 
shield materials. Of course the statistical varia- 
tions would be less in problems with less overall 
attenuation. 

The effect of statistical differences on the 
shield optimization can be reduced by following 
more packets. However, the computer times start to 
get excessive if this is the only approach used. 

It would be more fruitful in terms of the routine 
application of the program to expend some effort 
towards altering the importance sampling. 

The FASTER-III program has the capability of 
calculating optimal importance parameters based on 
partial derivatives of the variance. This feature 
can be used in determining better importance sam- 
pling parameters for shield optimization problems. 
In fact, the overall program efficiency could be 
improved if this feature was utilized on a wide 
variety of problems with the results being used to 
improve the built-in importance sampling models and 
parameters. 
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TABLE 1 


SPHERICAL REACTOR-SHIELD CONFIGURATION 


K levant 


COMPOSITIONS (10 2 " atome/cm 3 ) 

CORE 

SEPUICTO 

^38 

SHIELD 

MIX I 
SHIELD 

MIX 2 
SHIELD 

MIX 3 
SHIELD 

HoO+B 

SHIELD 

H 

0.01976 

0.0 

0.0 

0.0451 

0.0516 

0.0580 

0.0645 

Bo 9 

0.0 

0.120 

0.0 

0.0 

0.0 

0 0 

0.0337 

B 

0.0 

0.0 

0.0 

0.000671 

0.000766 

0.000862 

0.000958 

0 

0.01184 

0.0 i 

0.0 

0.0226 

0.02S8 

0.0290 

0.0 

AL 

0.0512 

0.0 

0.0 

1 0.0 

0.0 

0.0 

0.0 

Zr 

0.01744 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

u 2 * 

O.OOO979 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0*38 

0. 000078 

0.0 

0.0482 

0.01446 

0.00964 

0.00482 

0.0 


TABLE 3 

RESULTS OF FASTFE-III SHIELD OPTIMIZATION 
(0.25 mr/hr at 30 feet) 


TABLE 2 


RESULTS OF FASTES-III BASE LIME CALCULATIONS OF REACTCR 
SHIELD CONFIGURATIONS AT 30 FEET FROM CORE CENTER 


DOSE RATE CONTRIBUTIONS AT 30 FEET 
FROM COTIE CENTER 



375 ^ 

600 MW 

DOSE RATE 

REACTOR 

REACTOR 

COMPONENT 

(mr/hr ) 

.(»nr/hr) 


Photon Source Region 
Core 

0.009 

0.004 

Reflector 

3.5xlO' 6 

6.3xlO -6 

Depleted Uranium 

3.2xlo~ 5 

1.3xlO -5 

Mix 1 Shield 

0.018 

0.026 

Mix 2 Shield 

0.062 

0.075 

Mix 3 Shield 

0,017 

0.063 

Berated Water Shield 

O.Oll 

0.022 


Total Photons 

0.120+0.034 

0.187+0.054 

Neutrons 

0.020+0.002 

0.027+0.003 

Total 

0. l4o 

0.214 


Quantity 

37 

RE 

'5 MW 

ACTOR 

600 MW 
REACTOR 

Initial 

Final 

Initial 

Pina l 

Dose Rate (mr/hr) 





Photon 

0.120 

0. L26 

0. 187 

0.153 

Neutron 

0.020 

0. 124 

0.027 j 

0.097 

Total 

0. 140 

0.250 

0.214 

0.250 

Shield Weight (10 3 kg) 





Depleted U 

10.2 

12.6 

13.8 

0.0 

Mix 1 

71.2 

0.0 

89.2 

6.6 

Mix 2 

22.1 

52.4 

26.4 

52.4 

Mix 3 

16.1 

12.2 

19.0 

63.1 

Water 

§ Ll 

80.3 

97.7 


Total 

206.3 

157.5 

246.1 

j 207.4 

Shield Thickness (cm) 





Depleted U 

5.0 

6.1 

5-0 

0.0 

Mix l 

57.0 

0.0 

57.0 

7.0 

Mix 2 

15.0 

57.3 

15.0 

48.4 

Mix 3 

15.0 

13.5 

15.0 

51.4 

Water 

117.0 

120.8 

U7.0 

98.4 
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